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Abstract 



Branch points in the complex plane are responsible for avoided level crossings 
^ ■ in closed and open quantum systems. They create not only an exchange of 

the wave functions but also a mixing of the states of a quantum system at 
high level density. The influence of branch points in the complex plane on the 
low-lying states of the system is small. 
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I. INTRODUCTION 

The statistical properties of the shell-model states of nuclei around ^'^Mg are studied a 
few years ago [0] by using two-body forces which are obtained by fitting the low-lying states 
of different nuclei of the 2s — Id shell. In the center of the spectra, the generic properties 
are well expressed in spite of the two-body character of the forces used in the calculations. 
It is therefore not surprisingly that embedded random matrix ensembles are relevant for a 
description of eigenfunctions and transition strengths in many-body quantum systems |0-§1 . 
Some regularity in the ground state region of the spectra is obtained in recent calculations 
with the two-body random ensemble [§-0] which could be explained in the cases considered 

In the calculations with the two-body random ensemble, the two-body character of the 
forces in the Hamiltonian and the Pauli principle are taken into account: states with given 
total angular momentum J consist of A^ identical Fermi-particles distributed among a set of 
single-particle states, and interact via two-body matrix elements. By standard shell-model 
(SM) techniques, the Hamiltonian matrix is calculated in terms of the independent two- 
body matrix elements which are taken at random from a given distribution. The number of 
independent two-body matrix elements is small as compared to the total number of matrix 
elements. Due to its nearness to the SM, the two-body random ensemble is expected to 
be more appropriate than the Gaussian orthogonal ensemble for a study of nuclear spectra 
|I^,|TT[]. The level density calculated with a two-body random ensemble is, indeed, more 



realistic than that of the Gaussian orthogonal ensemble. It is nearly Gaussian in contrast 
with the Wigner's semi-circle law, that holds for the Gaussian orthogonal ensemble in the 
limit of a large number of states. 

In Ref. 0, the spectrum of A^ fermions distributed over M orbitals and the structure 
of the eigenstates is calculated with a two-body random ensemble by using the methods 
of scar theory. At the edges of the spectrum, significant deviations from random matrix 
theory expectations are found. The deviations are much smaller when a spin is added to 
the Hamiltonian (as in nuclear SM calculations). Due to the spin, the single-particle states 
split and the many-particle SM states may cross. The crossing of any interacting discrete 
states is, however, not free but avoided. Thus, the results of the calculations [0 with and 
without spin differ by the number of avoided crossings of their levels. 

It is the aim of the present paper to study the relation between avoided crossings of levels 
and the mixing of their wave functions in detail. The avoided crossings of the many-particle 
levels are a property of the whole system consisting of many particles. The correlations 
induced by them are a second-order effect (via the continuum) as compared to those induced 
directly by the two-body forces. Nevertheless, they may play an important role, especially 
at high level density. 

In section 2, the relation between branch points and avoided level crossings is discussed 
by means of a simple example. The correlations induced by the branch points in the complex 
plane are considered in section 3. Numerical examples are shown and discussed in section 4 
and the results obtained are summarized in the last section. 



II. BRANCH POINTS IN THE COMPLEX PLANE AND AVOIDED LEVEL 

CROSSINGS 

To every avoided crossing of discrete states, the corresponding crossing point in the 
complex plane can be found. This can be seen by considering the matrix 

/ei(a) V \ _i_ ( ci{a) \ , . 

[ V e2{a) J 2[ c^ia) ) ^'^ 

whose eigenvalues are 

^^, = E., - ^ Q, = '-^ ± ^(^1 - e.y + Av^ (2) 

with i, j = ±. Further, e^ = e^ — i/2 c^ {k = 1, 2) is the energy e^ of the state k continued by 
1/2 Cfc into the complex plane. The states with energies e^ are assumed to be unperturbed 
by any interaction with other states. Their energies may be traced as a function of a certain 
parameter a. Diagonalizing Ti gives the energies Sij of the eigenstates of 7i in which the 
interaction v between the two basic states 1 and 2 is taken into account. 

The eigenvalue equation (0) shows the following. When the non-diagonal matrix elements 
V (which describe the interaction between the two states k = 1,2) are real and different from 
zero, the two states i,j = ± can cross, as a function of the parameter a, only in the complex 
plane. The crossing point a = a™ is defined by £^+(a™) = S-{a'^^) and 

F'^' = {ei{a'') - £2(0")}^ + 4t;2 = . (3) 

This crossing point is a branch point in the complex plane according to Eq. (0). It lies at 

The non-trivial mathematical properties of these branch points in the complex plane as 
well as their relation to the higher-order poles of the 5* matrix are known for a long time. 



e.g. |T2|. The eigenfunctions $ij oiTC are bi-orthogonal. They are normalized |jl3|] according 
to 

($k)2 = (I>^,|I>^) = <5w (4) 

and therefore 

\^r\'' = {^r\^r)>1; \{^R'^R\^R)\ = \lm{^R,^n\^R)\>0. (5) 

Approaching a double pole, 

•IJ'P ^ ± 2 l>^P (6) 



see 



T4[, and 

bPia , ^^ . I /-T,bp I ;bp 



I'^S'r-oo; lW^«|$^^)|-oo, (7) 

see [|^. At an avoided crossing in the complex plane and in its neighborhood, the difference 
between |$_rP and ($i?)^ is large but finite. Numerical examples are given in [0,|l^]. The 



S matrix contains (^Ijl'l'/j) but not ($^^1$^^), i.e. it behaves smoothly also when the bi- 
orthogonality is large. This holds also at a double pole of the S matrix where two resonance 
states cross |T5l . 

The influence of avoided level crossings on the dynamics of a system under adiabatic 
conditions and on the correlations between the different states has been studied in atoms 



1^ , P^JT7| and open quantum billiards |[T8|JT9|1 . In these systems, it is possible to trace the 



turn-over of the avoided crossings into a real crossing (corresponding to a double pole of 
the S matrix) by tuning the system parameters. All the results obtained show that branch 
points in the complex plane may have an important influence on the spectroscopic properties 
of a quantum system although their number is of measure zero. 

III. CORRELATIONS INDUCED BY BRANCH POINTS IN THE COMPLEX 

PLANE 

The correlations between states are expressed usually by the mixing of their wave func- 
tions 

^n = Y.bR,<^' (8) 

i 

where the {^^} are the eigenfunctions of the Hamilton operator Tio of the unperturbed 
system. In the SM calculations, the functions {$°} are the Slater determinants. They do 
not contain any residual forces and represent therefore a natural basis set for the SM wave 
functions $^^. In the case of an open quantum system, the functions {$°} are the wave 
functions of the corresponding closed system. In the continuum shell model, these are the 
SM wave functions {$jj^}. In the case of the matrix (^, Tio is defined by f = and the 
{$°} are the eigenfunctions of Ho- 

The Schrodinger equation of an open quantum system with the Hamiltonian Ti can be 
written as 

{no-Sn)\^R) = -W\^n) 

= -Y,{^r,\W\^r)Y.{^r>\^r.)\^r.) (9) 

R' R" 

where the operator W describes the interaction of the two states R and R' via the environ- 
ment. In the case of the continuum shell model, this is the interaction via the continuum 
which is complex, as a rule [^,^. The source term (rhs. of Eq. (^) contains the {^r\^r/). 



Thus, the mixing of the wave functions via the continuum is related to the bi-orthogonality 
of the eigenfunctions of 7i, Eq. (H). Numerical examples for the bi-orthogonality, given in 
e.g. Refs. ||13| , p!5| , p!8| , p2| , show that the standard relations ($^1$^/) = Srr' are not fulfilled 
neither at the branch point in the complex plane nor at an avoided level crossing, see Eq. 
(^). Thus the source term of Eq. (^) contains terms being non-linear in the wave functions 
$/? when the quantum system is open and the ^r are complex. 

SM calculations are performed with an Hermitian Hamilton operator TC^^. The eigen- 
values and eigenfunctions are real. It seems therefore, that the source term in Eq. (j^) can 
not play any role in the SM calculations. This is, however, not true as will be shown in 



the following. The reason for this unexpected behaviour is, above all, the analyticity of the 
wave functions. The eigenvalues lying at negative energy show, as a function of tuning the 
system parameters, the same motion in the complex plane as those at positive energy. 

In the following, the motion of the eigenvalues by varying the system parameters will be 
considered in a language more fitting into the concept of the SM. As a matter of fact, SM 
calculations provide not only the eigenvalues Ei of the SM states but also the eigenfunctions 
<l>j. Knowing the SM wave functions ^f of the system consisting of A particles and those of 
the system with one particle less, the amplitudes ($f |$;^~^) of the spectroscopic factors can 
be calculated. Here, $^~^ is the wave function of a well-defined state of the residual system 
having one bound particle less than the initial system. These amplitudes depend on the 
quantum numbers of the separated particle. By multiplying the spectroscopic factor with 
an energy dependent common penetration factor, one gets the probability for the transition 
of the initial state i (with A bound particles) to the decay channel c (with one particle in 
relative motion to a certain state of the system with A — 1 bound particles). The separated 
particle may either be transferred (at negative energy) or emitted (at positive energy). In 
the last case, the product of spectroscopic factor and penetration factor is identified with 
the partial decay width (7^)^ of the state i into the channel c. Further, Fj = Y^dltY is the 
width of an isolated state. Thus, the most characteristic part of the width Fj of an isolated 
state is involved in the spectroscopic factor calculated from SM wave functions. 

As a result, the SM eigenvalues and eigenfunctions of the two bound systems with A 
and A — 1 particles, respectively, contain all information on the motion of the complex 
values Ei — ^Gt (where Gi is proportional to the spectroscopic factor) by varying the system 
parameters, although both Hamiltonians Ti.Q and Tig ~^ are Hermitian and their eigenvalues 
are real. The accompanying motion of the energies Ei is considered in many studies and 
avoided level crossings have been observed. The accompanying motion of the spectroscopic 
factors, however, is not considered up to now. 

IV. NUMERICAL RESULTS 

In order to illustrate the motion of the energies and spectroscopic factors as a function 
of a certain parameter, let us consider the matrix 7i, Eq. (|l]), with ei = 1 — a/2; 62 = 
a; Ci = C2 = and v real. By means of the parameter a, crossing or avoided crossing 
of the two states can be traced under different conditions, e.g. for different values of v. 
The results show crossing for v = and avoided crossing for all f 7^ at the critical 
value a'^^ = 2/3. The influence of v on the eigenvalues is small as long as v is small: the 
typical avoided-crossing behaviour can be seen (Landau-Zener effect, figure top left). The 
two states are exchanged, what is surely without any interest for statistical considerations. 
For larger v, the avoided crossing of the eigenvalues as a function of a is difficult to see 
in the trajectories: the eigenvalues are weakly dependent on a. The comparison with the 
eigenvalues for v = shows, however, the changes of the trajectories in the neighborhood of 
a^"^. The wave functions are exchanged also in this case at a = a^"^. The wave functions are 
mixed strongly near a'^^ for all v as can be seen from the right-hand side of figure |l|. Here 
the coefficients bfj of the representation $j = J^bij^^, Eq- (|), for v = 0.05, 0.5 and 1 are 
shown ($^ are the wave functions at f = 0). 



While the wave functions are mixed only in a very small region around the critical value 
a" for small f , they are mixed also at large distances \a — a'^^\ when v is larger (figure p. 
At high level density, avoided level crossings with other states appear, as a rule, before 
hi^i -^ 1 , &ijyi — i> is reached. As a consequence, the hi j contain components from more 
than just the two crossing levels. 

For illustration, the energies and wave functions of four states are calculated from 



n 



(4) 



/ei(a) V V "^ \ 

V 62(0) V V 

V V 63(0) V 
\ V V V 64(0) ) 



(10) 



and shown in figure ^ for the two values v = 0.005 and 0.03 (ci = 1 — a/3; 62 = 1 — 
(5/12) a; 63 = 1 — a/2; 64 = a). As long as v is small, the influence of the branch points 
consists mainly in an exchange of the corresponding two states and of the mixing of just 
these two wave functions in a region |a — a'^''| which is well separated from other analogous 
regions around other critical values. At larger v, however, the states are exchanged and 
mixed with all the other ones since the different regions |a — a^''|, corresponding to different 
branching points, overlap. 

In the special case considered, the regions |a — a^'"! in which fej j — > 1 , feij^i — *> is 
reached, are well separated for the different critical values a'jl^ as long as f < 0.02. The 
overlap of the different regions \a — a"\ starts at f ~ 0.02. Here, the mixing of the wave 
function of every state with those of all the other ones starts. It changes smoothly with 
increasing v > 0.03. Characteristic of this regime is that an unequivocal identification of 
the i and j in the mixing coefficients 6? ■ is impossible. In this manner, correlations between 
the wave functions of all the states may be created by the avoided crossings, more exactly, 
by the branch points in the complex plane (expressed mathematically by the source term in 
the Schrodinger equation (D). 

Due to the mixing, the wave functions of the states become more similar and the differ- 
ences of their spectroscopic factors in relation to a certain channel become less. That means, 
the level repulsion along the real axis is accompanied by an adjustment of the spectroscopic 
factors. This result coincides fully with the results obtained from a study of the motion of 
the poles of the 5" matrix in the neighborhood of an avoided crossing in the complex plane. 
Here, level repulsion along the real axis is accompanied by an adjustment of the widths 



14,18 



The results show that the energies of the states are almost independent of a and v, ii v 
is sufficiently large or (and) the level density is high. That means, the system is stabilized 
against (real) distortions of any type. This is true also for the distortions introduced by the 
real part of the coupling to the continuum (principal value integral) which does not vanish, 



generally [^. In any case, the states have lost their individual SM properties by spreading 
their spectroscopic features. This spreading is caused by the source term of Eq. (^) which 
continues analytically into the function space of the discrete states. It does no longer show 
features characteristic of the special type of the two-body forces used in the calculation. 

It follows immediately, that the states in the centre of a SM spectrum show features 
caused by the branch points in the complex plane, i.e. by the non-linear terms in the 
Schrodinger equation. The states at the border of the spectrum are less influenced by 

6 



avoided crossings with other states. This is true especially for the lowest states which may 
shift downwards by crossing only a few (or no) other states. 

V. SUMMARY 

In this paper, the influence of branch points in the complex plane onto the correlations 
between the states of a quantum system is studied. As a quantitative measure of the corre- 
lations, the degree of mixing of the wave functions, i.e. the values b^i in the representation 
(I) is used. The results are the following. 

- The exchange of two discrete levels which avoid crossing can be traced back to the 
existence of a branch point in the complex plane. Here, the levels have equal energies 
and widths, their wave functions are linearly dependent, and the Schrodinger equation 
contains non-linear terms. 

- The correlations induced by the branch points in the complex plane appear between 
the discrete states of a closed system in the same manner as between the resonance 
states of an open system. 

- The number of avoided level crossings at high level density is much larger than usually 
assumed. Most important correlations are introduced even by those avoided crossings 
which are difficult to identify at first sight. 

It follows further from these results that the correlations induced by the branch points 
in the complex plane are the larger the higher the level density is. At the edges of the 
spectrum, branch points have only a small infiuence on the spectrum. The properties of 
these states are therefore determined, at least partly, by two-body forces. 
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FIG. 1. The energies Ei (left) and mixing coefficients bf , (right) for two states i = 1,2 as a 
function of a for v = 0.05 (top), 0.5 (middle) and 1.0 (bottom). The circles (i = 1) and stars 
{i = 2) on the left-hand side show the energies of the two states at f 7^ and the full lines show 
their energies for v = (which cross at a = a" = 2/3). The full lines on the right-hand side are 
the coefficients bfj^^ while the dashed lines are the bf ■,^. Note the different scales (top, middle 
and bottom) of the Ei and 6?- axes. 
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FIG. 2. The energies Ei (left) and mixing coefficients 6?- (right) for four states as a function 
of a at f = 0.005 (top) and 0.03 (bottom) for ah couphng matrix elements. The full lines on the 
left-hand side show the energies of the states at f 7^ and the dash-dotted lines show their energies 
at ?; = 0. It is bfj^^ > hf -^^ ( ^f ,=j = bf ,^^ = 0.5 at o") for all i at small v (top right) as in figure 
||. At large v, every state is mixed with all the other ones and an unequivocal identification of the 
i and j in the mixing coefficients 6? • is impossible. 
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